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Abstract 

We present a numerical technique employing the density of partition function zeroes (i) to 
distinguish between phase transitions of first and higher order, (ii) to examine the crossover 
between such phase transitions and (iii) to measure the strength of first and second order phase 
transitions in the form of latent heat and critical exponents. These techniques are demonstrated 
in applications to a number of models for which zeroes are available. 
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1 Introduction 



A long standing concern of statistical physics is how best to distinguish between phase transitions 
of first and second order (for recent reviews, see [|], First order phase transitions involve 

the coexistence of two distinct phases and are characterised in the infinite-volume limit by the 
existence of a discontinuity in the first derivative of the free energy. For temperature (energy) 
driven first order phase transitions this discontinuity in the internal energy is the latent heat while 
for field (order parameter) driven first order phase transitions the discontinuity is the spontaneous 
magnetization. The specific heat and the magnetic susceptibility - the second derivatives of the 
free energy with respect to the temperature and external field respectively - then exhibit delta 
function singularities in these cases. This contrasts with the case of a phase transition of second 
order, where the appropriate second derivative of free energy diverges along with the correlation 
length. This divergence is characterised by critical exponents labelling the universality class of the 
system. The infinite volume system does not anticipate the onset of a phase transition of first 
order as the transition region is approached. There is therefore no critical region and no critical 
exponents in the usual sense of second and higher order phase transitions. 

Finite volume systems do however anticipate the presence of the phase transition as the thermo- 
dynamic limit is approached for both of the above scenarios. Numerical methods of identifying 
the order of a phase transition exploit this anticipation by including finite size scaling (FSS) of 
extrema of thermodynamic quantities which are singular in the thermodynamic limit at the tran- 
sition point. The counterparts of these singularities in finite systems are smooth peaks, the shapes 
of which depend on the order and the strength of the phase transition. The theory and use of FSS 
is by now well established for both first @, |, |, @, § and second order g |DJ scenarios. 



An increasingly popular alternative to the analysis of these extrema is the use of the FSS behaviour 



of partition function zeroes |12| , 13, 14, [Lj, |16[] . For field driven phase transitions one is interested in 
the Lee- Yang zeroes in the plane of complex external magnetic field [|l^]. For temperature driven 
phase transitions the Fisher zeroes in the complex temperature plane are relevant [13]. In this 



case, FSS of the first Fisher zero yields the correlation length critical exponent v. By formally 
identifying v with 1/d in the first order case, FSS of this first zero can often be used to discriminate 
between the two types of phase transition. The strength of the transition, on the other hand, has 
to our knowledge heretofore not been determined from properties of the zeroes, but only by direct 
measurement of the latent heat or the interface tension. Analogous statements apply in the field 
driven case. 

Here we wish to present an alternative method to investigate the strength of phase transitions. 



This method involves the density of zeroes - a quantity directly proportional to the latent heat in 
the first order case. The emphasis in this paper is on interpretation rather than new simulations 
and our motivation comes from recent experience with numerical studies of the density of zeroes 
[17, IS | coupled to the ongoing discussions regarding first order phase transitions 0. The rest of 
the paper is organised as follows. The properties of partition function zeroes are discussed in Sec. g. 
In Sec. |H we present the results of our analysis applied to previously published data, and Sec. || 
contains our conclusions. 



2 Partition Function Zeroes 

The study of partition function zeroes was pioneered by Lee and Yang [^] in the complex external 
ordering variable for field driven phase transitions. Study of zeroes in the complex temperature 
plane was initiated by Fisher ifDf l. For convenience, one often refers to the zeroes in the complex 
external field strength plane as Lee- Yang zeroes and those in the complex temperature plane (in 
the absence of an external field) as Fisher zeroes. The emphasis in this paper is on Fisher zeroes, 
although the techniques can also be applied to Lee- Yang zeroes. 

Itzykson, Pearson and Zuber [|l4j] gave the FSS of partition function zeroes for <i-dimensional systems 
below the upper critical dimension. They suggested a compact description of scaling which for the 
j th Lee- Yang and Fisher zero respectively reduces to (for large j) 

d+2-T) _ J_ 

Here, L is the system size, rj is the anomalous dimension, t = T/T c — 1 is the reduced temperature 
which is zero in the first formula, and h denotes the external field which is zero in the second 
formula. The integer index j increases with distance from the critical point. The FSS at and above 
the critical dimension was determined in and [15] respectively. This behaviour of partition 



function zeroes offers a (generally very accurate) way to determine (i) the critical temperature and 
(ii) the critical exponents which describe the thermodynamic critical properties of the system. 

Previous numerical studies of partition function zeroes mostly employed the following strategies: 

(i) Determination of the distribution of zeroes to study the complex parameter phase diagram and 
to check (in a qualitative way) if the zeroes appear to cross the real axis - evidence for the existence 
of a phase transition (l|, |o|, ||, ||, ||, ||, ||. 

(ii) Quantitative measurement of (finite size) scaling behaviour of those zeroes with the smallest 
imaginary parts to determine the critical exponents [16, 26, 27, 28, p^, 3C, 31, |35, 33, 34 with 



possibly similar analyses applied separately to higher index zeroes |3l| or determination of impact 



angle ||36|, |37l, |3 



(Hi) Fixed volume plots for Lee- Yang zeroes indicating onset of non-zero density of zeroes at first 



order phase transitions |24|, |39[ . 



Despite an early study [4C], it was long considered difficult if not impossible to estimate the density 
of zeroes from their distribution for a finite lattice plf|. I n recent years, however, there have 



been some attempts to extract the density of zeroes from numerical studies [17, 41 1. In view of 
the increasing importance attached to this approach, we wish to suggest an appropriate way this 
should be done. 

To this end, we present a novel and powerful technique with which (i) the order of the phase 
transition and (ii) the strength of the phase transition can be determined. This means determination 
of the latent heat in the first order case and measurement of the critical exponent a of the specific 
heat in the second order case. The method also elucidates crossover between second and first order 
transitions. We study the cumulative distribution of zeroes - curves involving L and j collapse - 
from which the density of zeroes and hence the order and strength of the phase transition can be 
determined. 



2.1 Density of Zeroes 

We start with the factorised representation of the partition function, 

Z L (z) = A(z)H(z-z j (L)) , (2.2) 
j 

where z stands generically for an appropriate function of temperature (in the Fisher case) or field 
(in the Lee- Yang case), L is the linear extent of the lattice and A(z) is a smooth function that 
never vanishes. The free energy density follows as 

f L {z) = ^ d \nZ L {z) = ^ d \nA{z) + ±- d Y J ^{z-z j (L)) . (2.3) 

j 

The first term on the right contributes only to the regular part of thermodynamic functions and 
we henceforth drop it. The remainder, which we refer to as gives rise to singular behaviour. 

Following Abe f42|] , we assume the zeroes, Zj, (or at least those close to the real axis and hence 
determining critical behaviour) are on a singular line for large enough L, impacting on to the 
real axis at an angle (p at the critical point z = z c . The singular line is parameterised by z = 
z c + r exp (i(p). Define the density of zeroes as (with Zj = z c + rj exp (i(f)) 

g L (r)=L- d ^(r-r,{L)) . (2.4) 
j 



The free energy is 

fl(z) = [ R g L (r)ln (z-z c -re llfi )dr + c.c. , (2.5) 
J o 

where c.c. means complex conjugate and R is some appropriate cutoff. The cumulative distribution 
function of zeroes is defined as 

G L (r) = f g L (s)ds . (2.6) 
Jo 

For a finite size system, therefore, this is a step function with 

G L {r) = j- d if re( rj ,r j+1 ) . (2.7) 



It is natural to assume that at a zero, this cumulative density is given by the average [12, 15, 43] 



Gl ^ ) = \j- ■ (2 - 



In the thermodynamic limit and for a phase transition of first order Lee and Yang |12j already 
showed that the density of zeroes has to be non-zero crossing the real axis, 

9oo(r) = goo(0) + ar w + . . . . (2.9) 

This corresponds to the cumulative distribution of zeroes 

G OD (r)=g OD (0)r + br w+1 + ... , (2.10) 

where the slope at the origin is related to the latent heat (or magnetization) via |T^] 

Aeoc 5oo (0) . (2.11) 



Abe [42 1 has shown that the necessary and sufficient condition for the specific heat at a second 



order phase transition to have the leading critical behaviour C ~ t a , is (see also |44j| ) 

5oo (r) oc r 1 '" . (2.12) 
If a = 0, as is the case in the d = 2 Ising model, ( |2.12| ) leads to a logarithmic divergence in the 



specific heat [42]. The corresponding expression for the cumulative distribution (integrated density) 
is 

Goo(r)ocr 2 - a . (2.13) 

Thus while the scaling behaviour of the position of the first (few) zeroes in the complex temperature 
plane can be used to identify the density of zeroes gives the strength of the transition. A plot 
of Gx(rj) = (2j — l)/2L d against rj(L) should (i) go through the origin, (ii) display L- and j— 
collapse and (Hi) reveal the order and strength of the phase transition by its slope near the origin. 



2.2 Finite Size Scaling 



The above density approach yields new insights into traditional finite size scaling which emerges 
quite naturally from it for both first and second order phase transitions. Traditionally, FSS is based 
on a hypothesis whereby the only relevant length scale of the system is the correlation length. The 
FSS hypothesis can, in turn, be justified on renormalization group grounds. 

Alternatively, and in the second order scenario, one may start from (|2.8| ) and ( [2.13 ) with fixed j 



and identify Gl and Goo for large enough L. The distance of a zero from the critical point is 

'2j_i\V(2-«) _ /2j-l\ 1 /" d 



"<®~\nr) =\rir) (2 - 14) 

from hyperscaling. This is the standard formula for FSS of Fisher zeroes at a second order phase 



transition implied by [14]. It also recovers the compact expression (|2.1|) for high index zeroes. 



From this, FSS for all thermodynamic quantities can be found fllq , 42, ^4[. For the specific heat, 
for example, (|2.3j ) yields 



CUz) = ^TT^y ■ ( 2 ^) 



j 

At z = z c , and assuming that pseudo critical behaviour is dominated by the lowest few zeroes [ pf , 
one finds the FSS of specific heat is 

C L (z c )~L- d (z c -z j )- 2 ~L- d L 2 /" = L<*/» . (2.16) 

Thus the usual formulae for leading order FSS are recovered. 
Also, in the first order case, (|2.8| ) and fl2.10|) , with fixed j give 

rj {L)~L- d , (2.17) 

and hence Cl(z c ) ~ L d , explaining the formal identification of v with 1/d alluded to in the intro- 
duction. 



3 Analysis of Various Models 

In the following examples we perform fits to the cumulative density of zeroes allowing for first or 
second order behaviour, 

G(r) = ai r a2 + a 3 . (3.1) 

where we also allow for an additional parameter 03. The criterion for a good fit, apart from good 
data collapse (in L and j), is that 03 be compatible with zero. In fact, a value of 03 inconsistent 



with zero indicates the absence of a phase transition, for, if 03 > 0, the zeroes have already crossed 
the real axis corresponding to the broken phase and if 03 < the zeroes remain away from the real 
axis which is the situation in the symmetric phase. A first order phase transition is then hinted 
at if a% ~ 1 for small r (constant slope near the origin). Then the latent heat is proportional to 
the measured slope a\. A value of 02 larger than (and not compatible with) 1 is indicative of a 
second order phase transition. In this case, the strength of the transition is given by a = 2 — 0,2- 
The cumulative distribution of zeroes is given by ( |2.8| ) and the parameter r may be taken to be 
the imaginary part of the position of the j th zero. That is, in practice, we fit to 

^l = G L (r j ) = a 1 (lmz j (L)r+a 3 . (3.2) 

A straightforward error analysis is not possible because of the semi-discrete nature of this equation. 
We employ a procedure adapted from and which requires some explanation. 

We observe that (f2.8j) arose from the natural assumption that the density associated with a zero 



in the finite volume case is obtained from averaging the exact formula (2.7). If one relaxes this 
assumption, and, since the cumulative density is monotonic, its actual value cannot deviate from 
the mean position fl2.8| ) by more than ±l/2L d . Denote the data point corresponding to the j 
zero of the size L lattice by G° bs (L). Inspired by |45], initially, assign an error (Jj{L) = cx ar b/£ d to 
this data point, where cr ar b is a fixed arbitrary number. Set 

X? = E — S(L) 7nf P(L))2 = E ^(G?°(L) - Gf*{L)f , (3.3) 

L,j ' L,j °arb 

where G^ xp (L) is the expected density value coming from the model ([O]) and the sum runs over 
the lattice sizes and indices used in the fit. One proceeds in the usual way to fit the parameters 
in ( |3.lD by minimizing Xi- The errors associated with a,i resulting from this fit are denoted daf rh . 

Assume, now, each data point has, in fact, an error a/L d . The corresponding chi-squared is 

Xl = E ^(Gf\L) - Gf*{L)f = a -^ X \ . (3.4) 

L,j 

Assuming, the model fits well, this should be close to the number of degrees of freedom A^of- In 
this case, the error assigned to each point becomes 

a 2 = al ThX i/N do{ = xi/N doi , (3.5) 

if cr ar b is chosen to be unity. Furthermore, the true errors associates with the fit parameters are 

d ai = —daf h = adaf h (3.6) 

Carb 

if farb = 1. One notes, however, that as in [45], this approach prohibits an independent goodness- 
of-fit test. 
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Figure 1: Distribution of Fisher zeroes for the d = 2, q = 10 Potts model which has a strong first 
order phase transition. The symbols x, +,A,<>,n, and O, correspond to j = 1,2,3,4,5, and 6, 
respectively, obtained on lattices of size L = 16-64. 

3.1 The Potts Model 



The g-state Potts model is the classic testing ground for the analytical and numerical study of first 
and second order phase transitions. Despite the lack of a complete solution for general q, some 
exact results are known. In two dimensions, the model has a second order phase transition for 
q < 4 and a first order phase transition for q > 4. In terms of the inverse temperature, (3 = 1/ksT, 
the critical temperature is given by (3 C = ln(l + yjq) and the exact values of the latent heat and a 
number of other thermodynamic quantities are known p6(| . In the three dimensional Potts model, 
it is numerically established that there is a second order phase transition in the Ising case (q = 2) 
and a first order phase transition for q > 3. For a review on the Potts model, see [47]. 

In the following examples we use the zeroes published in a number of papers. If, for the Fisher 
zeroes, the parameter z used is exp(— k(3) where k is a trivial scale parameter depending on the 
model under consideration, then the latent heat (for a first order phase transition) is 



Ae = Kz c 2-Kg(0) , 



(3.7) 



where z c = exp(— k(3 c ) is the coupling at the phase transition point and g(0) is the slope of the 
cumulative density plots. 



Table 1: Fits of the d = 2, q = 10 Fisher zeroes to the cumulative distribution G = g(0) x Imz. 
Here, the N lowest zeroes from the L = 16-64, j = 1-4 data are used. As the origin is approached, 
the resulting estimates of the latent heat approach the exact value 0.6961. 



N 


24 


20 


16 


12 


08 


04 


5(0) 


0.501(8) 


0.486(4) 


0.479(3) 


0.471(2) 


0.469(2) 


0.463(1) 


Ae 


0.756(11) 


0.734(6) 


0.723(4) 


0.711(3) 


0.708(2) 


0.698(2) 



The d = 2, q = 10 Potts Model: Villanova |3C] has listed the first six Fisher zeroes in the 
z = exp (— 0) plane (i.e., k = 1) for the two-dimensional 10-state Potts model for lattice sizes 
L = 4,8, 16, 24, 32, 38, 48, and 64 (see also [pl| ). In |3l| a conventional FSS analysis was applied 
to the first index zeroes only. The leading scaling behaviour using the three largest lattices yielded 
an acceptable fit and v = 2.0026(11). Including the leading additive correction term gives an 
acceptable fit for L = 24, 32, 38, 48, and 64 with the result v~ l = 2.0028(17). This is good evidence 
for v = 1/d and hence a first order phase transition. However, the determination of the lattice size 
above which FSS (or FSS with corrections) sets in is, by necessity, somewhat arbitrary. Indeed, 
when one extends the analysis to higher index zeroes (FSS applied to each index separately), one 
finds that when corrections are ignored, no two-parameter fit gives an acceptable result. Allowing 
for corrections however (fitting to four parameters) may yield acceptance indicative of a first order 



phase transition. These higher index zeroes were not however analysed in [30, 31]. 

Our analysis of the density of zeroes begins with Figure [l] where all six zeroes are plotted for 
L = 16-64. We do not plot the results from smaller lattices simply for clarity - including them 
will only give some more data points to the right of the figure and will not affect the slope in the 
region of interest, namely near the origin. Figure [l] displays excellent L- and j- collapse indicating 



that ( ]2.8| ) is the correct functional form of the density of zeroes. Fitting (3.1) to the L = 16-64, 
j = 1-4 points gives a± = 0.84(2), 02 = 1.10(1) and 03 = 0.00004(1) strongly indicative of a first 
order phase transition. Fixing 03 = 0, and fitting for the two remaining parameters again gives 02 
close to 1. Indeed, such a fit to the lowest four data points gives 02 = 1.008(6). Fixing 02 = 1 
(and as = 0) and applying a single parameter fit to the full data set yields g(0) = a± = 0.501(8). 
Further fits close to the origin yield the slopes indicated in table [l]. From (|3.7| ), the latent heat is 
Ae = 27rg(0)/(l + \/l0) and these corresponding values are also listed in the table and appear to 
converge to the exact value 0.6961. 
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Figure 2: Distribution of the L = 18-36 Fisher zeroes for the d = 3, q = 3 Potts model which 
has a weak first order phase transition. The symbols O, □, and O correspond to j = 1,2, and 3, 
respectively. 

The d = 3, q = 3 Potts Model: It is generally accepted that the three dimensional three 



state Potts model has a first order phase transition, albeit a very weak one [48|. This is therefore 
a typical model used to test new methods to discriminate between phase transitions of first and 
second order. A list of the first few Fisher zeroes (in the z = exp (—3/3/2) plane) for the d = 3, q = 3 
Potts model for L = 10-36 can be found in fjO] , j3~lfl . In fact there is some ambiguity regarding the 
labelling of the higher zeroes in [ ]30| , |3l| as can be seen by comparing their distributions for different 
L. To be cautious, therefore, our analysis only employs the zeroes with index j = 1-3. Again in 



[30, 31], v is determined by applying FSS to the first index zeroes only. However, even with the 



first index zeroes, one finds rather bad fits for any combination of more than two lattice sizes and 



these are not much improved upon by allowing for corrections. The way out taken by [30, 31] is 
to measure v using the first index zeroes for pairs of consecutive lattice sizes in order to convince 
oneself that v approaches 1/d as L increases. Their result using only the two largest lattices is 
z/ _1 = 2.955(26). Thus the evidence presented points towards a first order phase transition, with, 
however, the caution that "with the presented trend from other lattices, one does not feel confident 
to exclude accidental agreement" . 

Our density analysis is presented in Figure |L A 3-parameter fit to all data yields = 0.000005(2) 
and becomes even closer to zero as the fit is restricted closer to the origin. Clearly the slope is non- 



zero near the origin - the signal of a first order phase transition. In fact, a 2-parameter fit to the 
data corresponding to L = 22,24,30,36, j = 1 yields a\ = 0.070(8), 02 = 1.06(2). Accepting that 
the plot is in fact linear near the origin, and fitting for the slope only gives g(0) = a± = 0.0454(9). 
Using f3 c = 0.3670 §0], ||] and Ae = (3/2) exp [-(3/2)/3 c ]2tt 5 (0), from Q with k = 3/2, we 
find that the corresponding latent heat is 0.247(5) comparing well with 0.2409(8) from [[}(], and 



with 0.2421(5) from the more sophisticated analysis of 48]. 



The d = 3 Ising Model: The Fisher zeroes in the z = exp(— 4/?) plane for the three dimensional 
Ising model with L = 4, j = 1-7 have been determined exactly in [pcf ] and are listed in [37] where 
the zeroes for L = 5, j = 1-4 are also determined. We also use the zeroes for L = 7, j = 1,2 
from ||, those for L = 6, 8, 10, 14, j = 1, 2, 3 which are listed in p0| , |32"| and j = 1 zeroes up to 
L = 32 recently obtained in [49]. The cumulative density of zeroes (using all of the above data) 
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Figure 3: Distribution of Fisher zeroes for the d = 3 Ising model which has a second order phase 
transition. The symbols x,+, □, and O correspond to j = 1,2,3, and 4, respectively and the 
lattice size ranges from 4 to 32. 



is plotted in Figure y. Fitting the full set of L = 4-32, j = 1-3 data to (|3.1[) indicates a second 
order phase transition with a<i = 1.81(3), 03 = —0.00001(1). Accepting the plot goes through the 
origin and applying a 2-parameter fit to the L > 6 data gives ai = 1.874(42) which corresponds to 
a = 0.126(42). Fitting closer to the origin (using the six data points corresponding to L = 10-32, 




Figure 4: Distribution of j = 1 Fisher zeroes for the Ising model in two dimensions which has a 
second order phase transition. The lattice sizes range from L = 4 to L = 64. 

Table 2: Fits of the d = 2, q = 2 Fisher zeroes to the cumulative distribution G = a\ x (Imz) a2 . 
As the origin is approached, the parameter a-i approaches 2. 



L 


4-64 


10-64 


16-64 


20-64 


24-64 


32-64 


36-64 


40-64 


48-64 


a-2 


1.924(9) 


1.958(4) 


1.968(3) 


1.972(2) 


1.975(2) 


1.978(1) 


1.980(1) 


1.981(1) 


1.982(1) 



j = 1 gives a>2 = 1.879(2) or a = 0.121(2), close to the expected value ~ 1/8. Note that the 
value of v was measured to be 0.612(22), 0.6295(10), and 0.6285(19) in references ||, 0, and 



[30, 32]. The corresponding a values from hyperscaling are a = 2 — vd = 0.164(66), 0.1115(30), 
and 0.1145(54), respectively. A weighted mean value of v from various literature estimates for the 
d = 3 Ising model is given in |5(| as 0.63005(18) and corresponds to a = 0.10985(54). 



The d = 2 Ising Model: The (numerically) exact Fisher zeroes for the Ising model in two 
dimensions for large lattices with periodic boundary conditions (up to L = 64) have only recently 



been determined in [51]. Their distribution with z = exp(— 4/3) is plotted in Figure | 



The results from two parameter fits to the L-ranges available are given in table Fitting near 
the origin gives ai = 1.982(1), close to 2 which corresponds to a = 0. This value for a is both 
sufficient and necessary to attain the logarithmic singularity of the specific heat in the d = 2 Ising 
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Figure 5: Distribution of Fisher zeroes for the d = 3, SU{2>) pure gauge model which has a first 
order phase transition. The insert is a blow-up of the region near the origin. 



model ||. 



Although only the first index zeroes are available in this case, the cumulative density method is still 
useful. The conventional FSS technique can only yield the ratio a/v from the analysis of specific 
heat. Furthermore, here the specific heat leading scaling behaviour is complicated by logarithms 
while the density and cumulative density are free from such terms [42|. Also, FSS applied to the 
zeroes themselves (in the conventional manner) yields the correlation length critical exponent v, 
which is related to a in cases where hyperscaling holds. The cumulative density method is thus 
the only method by which a can be measured directly from finite size data. 



3.2 Lattice Gauge Theories 



Since lattice gauge theories may be viewed as statistical systems (see, e.g., [52[ for a review), we 
complete this section by presenting some applications of our method to published zeroes of such 
theories. 



The d = 3, SU(3) Model: Alves, Berg and Sanielevici [ 33 1 have analysed the SU(3) finite 
temperature deconfining phase transition using the FSS of the lowest Fisher zeroes for Ltl? lattices 
(see also [^]). They determined the zeroes in (5 for Lt = 4 on lattices of size L = 4, 6, 8, 14, 16, 20, 
and 24. Their conventional FSS analysis gives the "illusion of a second order phase transition" when 
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Figure 6: Distribution of Fisher zeroes for the four dimensional abelian surface gauge model which 
has a second order (mean-field like) phase transition equivalent by duality to the four dimensional 
Ising model. The □ and O correspond to the j = 1 and j = 2 index zeroes respectively. 

lattice sizes L = 4, 6, and 8 are involved in the fits. Claiming FSS only sets in for larger lattices, 
their best result is v = 0.35(2) for L = 14-24. This value is compatible with v = 1/d = 0.33 and 
thus indicative of a first order phase transition. Figure ^ is a plot of the cumulative density of zeroes 
for all lattice sizes. The insert is the same plot using only L = 14-24. The figure, clearly supportive 
of a non-zero slope through the origin, justifies restriction of the analysis to these largest lattices 
and thereby elucidating the procedure of deciding where FSS sets in. This slope is 0.0121(3), which 
gives a latent heat of 0.0760(19). This is fully consistent with the result for the latent heat coming 
from conventional methods |^] which is 0.0758(14). 
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The d = 4, Abelian Surface Gauge Model: In the abelian surface gauge model, abelian 
gauge variables are assigned to plaquettes of the lattice. Models such as this are of relevance to 
supergravity theories amongst others. It is known that the d = 4 version is dual to the d = 4 Ising 
model, which, up to logarithmic corrections has mean field critical exponents p6| . One therefore 
expects the d = 4 abelian surface gauge model also to be characterised by mean field exponents 
(up to corrections) with a = 0,v = 1/2. 



This model was analysed in [53| where the first and second Fisher zeroes for lattices of size L = 3 
to L = 12 are listed. A conventional FSS analysis applied to the first index zero yields the best 



estimate of v = 0.469(17) from the two largest lattices (L = 9 and L = 12). It is remarked that 
inclusion of the smaller lattices worsens the fit and drives v away from 1/2. The appearance of a 
bimodal structure in the energy histograms was discussed in [|53[ as a spurious indication of a first 
order transition, at odds with the critical exponents which clearly favour second order. 



A fit of the data of 53] to fl3.1| ) yields 02 incompatible with unity (see Figure |6|). In fact, a two 
parameter fit near the origin (using the lowest 6 zeroes) yields 02 = 1.90(9) corresponding to 
a = 0.10(9), compatible with zero. Note, again, that only the region near the origin in Figure |6]is 
of interest. The slope at the origin is not compatible with a first order transition and the bimodal 
structure of the energy histograms observed in |^] can only be a (so far unexplained) finite size 
effect which merits further investigation. 



4 Conclusions 

The two basic questions to be addressed in the numerical study of phase transitions are (i) iden- 
tification of the order of the transition and (ii) determination of physical quantities such as the 
transition point and critical exponents or the latent heat. 

We have presented a method to determine the order and strength of phase transitions based on an 
analysis of the cumulative density of partition function zeroes. From the qualitative behaviour of 
this cumulative density we can distinguish between first and second order while from the quantita- 
tive details we can extract the latent heat in the case of first order transitions and the specific heat 
exponent a in the second order case. This should be compared with traditional finite size analyses 
where only the combination a/u can be extracted directly. 

Our method meets with a high degree of success even in the borderline cases of the d = 3, q = 3 Potts 
model and SU(3) finite temperature lattice gauge theory where the distinction between first and 
second order phase transitions is quite difficult. There, the conclusions based on the conventional 
FSS analysis are tenuous or ambiguous. Furthermore, the method leads to alternative insights into 
statistical physics methods in general and especially illuminates the origin of finite size scaling. 

Acknowledgements: R.K. would like to thank Wolfgang Grill for his hospitality during an 
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